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SUMMARY 


A theory is presented vhich appears capable of predicting the 
attenuation of hroad -hand sound in soft vailed (lov acoustical impe- 
dance) circular ducts. 

The axially symmetric separable solutions to the linear vave 
equation in clyindrical coordinates were obtained. It vas assumed 
that the steady velocities and pressure fluctuations vere small. The 
first ten characteristic function solutions vere fit to form a travel- 
ing plane pressure vave at the lined duct entrance by means of a 
Fourier-Bessel series expansion. The eigenvalues are complex n umb ers 
due to the finite vail impedance boundary condition. 

The present theory, in contrast to the approximate theory valid 
for hard vails, shovs a strong dependence of sound attenuation on 
frequency. This dependence occurs through the ratio of the duct 
diameter to sound vave length. 

At high frequencies focusing of the sound vave tovard the duct 
axis occurs. The acoustic energy can remain high, but at the outside 
vail, vhere the energy must be absorbed, the pressure oscillations may 
be small. 

The theory vas used to calculate the noise reduction obtained 
vith a perforated plate liner on a J-65 engine inlet. The liner de- 
sign vas based upon an approximate theory valid for very hard vails. 
The present theory shoved excellent agreement vith experimental 
measurements and vhy the design based on approximate theory failed to 
provide high noise reduction. 

SYMBOLS 

A area, ft^ 

Bj constant (see eq. (5) and (13)) 

C speed of sound, ft/sec 

0^ specific heat at constant pressure, Btu/lb m /°R 

D duct diameter, ft 

D/A duct diameter to vavelength ratio (see q belov) 

dB max maximum possible sound power attenuation for a given q 
and L/D, decibels 


E 


acoustic energy, ft lb^/sec 
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jk 


n 

L/D 

L 

P J 

P* 
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R j 


ft 


fc w 


T 

t 


X 

Z 

CXj 

J 

s 

n 


e j 

X 

p 

Po 


acoustic energy contribution of j 
velocity modes, ft lb^/sec 

sound frequency, Hz 

Bessel function of the first kind, 

duct length to diameter ratio 

duct length, ft 

j mode of acoustic pressure, psi 
complex conjugate of P, psi 


acoustic pressure. 



pressure 


of order 


and 


n 


k 


axial 


modulus of j complex eigenvalue 
gas constant, ft Ib.j./l'b /°R 
radial space coordinate, ft 
duct radius, ft 
temperature, °R 
time, sec 


acoustic velocity, 



ft/sec 


j mode of acoustic velocity, ft/sec 

axial space coordinate, ft 

liner acoustic impedance, Ib ffl /ft^/sec 

complex eigenvalue of j radial characteristic solution 

(aj - R j elej > 

dimensionless liner acoustical impedance (z/p C) 

Morse's frequency parameter, duct diameter in wavelengths 
r a>/jtC or D/A 

argument of j complex eigenvalue 
sound wavelength, ft 
density, lb m /ft 3 
average density, Xb m /ft^ 
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cTjjTj Morse's transmission parameters (ref. 2) 

0 dimensionless acoustical resistance, (Re (£) 

X dimensionless acoustical reactance, 4m U) 

circular frequency, rad/sec 
( )* complex conjugate 

1. INTRODUCTION 


The solution to the linearized acoustical wave equation in circu- 
lar ducts is well known and the method of applying boundary conditions 
to the characteristic solutions to represent a complicated waveform 
entering a soft walled duct is thoroughly established, (refs. 1 and 2). 
Due to computational difficulties, an approximate solution to the duct 
wall impedance equation for the first mode of vibration has been made 
in the past, (refs. 1, 3, and 4). The resulting solution is an ap- 
proximation to the transmission of a plane wave in a circular duct 
with a nearly hard wall. 

The purpose of this paper is to investigate the complete solution 
of the linearized wave equation with associated boundary conditions 
such that the solution is valid for ducts with very soft walls (i.e., 
acoustically lined ducts). The solution is specialized to that of a 
plane wave entering the duct but could be applied to any axisymmetric 
wave form. 

The theory was applied to broad-band or random noise in the fol- 
lowing way. At each frequency in the broad -band spectrum the random 
noise was assumed to enter the duct as a plane traveling pressure 
wave. The attenuation of this pressure wave was calculated at each 
frequency and compared to that obtained experimentally in the inlet of 
a J-65 turbojet engine. The comparison between experimental and 
theoretical sound power attenuation was excellent. 

2. THEORY 


Differential Equation and Solutions. The equations of motion 
for the air in the duct are: 


II 
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Viscous and heat transfer terms have been neglected. The air 
in the duct is assumed to behave as an ideal gas. If it is assumed 
that the pressure perturbations are small compared to the average 


pressure and the steady velocities are negligible, equation (l) re- 
duces to: 


- -po c2v • v 1 

p o| J 


Equations (2) can he combined to obtain the wave equation, 

p = ± ^2 

C 2 St 2 


(3) 


and if the time dependence in pressure is of the form e^ 0 ^ then 



(4) 


The solution of equation (3) in cylindrical coordinates with 
axial symmetry is: 

P j - it) expjtot-f (OJ + iTj)xj (5) 

The subscript j is used in equation (5) since there are an 
infinite number of solutions of this form. 

The radial eigenvalues must be determined from the bound- 

ary condition at the wall (r = r w ) and can be expressed as: 

_ i0n . v 

aj = Rj e “ (6) 

The constants B< njust be determined from the initial pressure 
wave at the duct entrance (X = 0). 

It should be noted that the argument of the Bessel function and 
thus the Bessel function itself is a complex number. The methods used 
for calculation of the complex Bessel function are presented in the 
appendix. 

Application of Boundary Conditions. Each solution (eq. (5.)) must 
satisfy the radial impedance condition at the outer wall: 



Equation (4) is used to obtain the radial velocity (v .). Equation (7) 
can then be expressed as : " 
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where 


r) = 



2r w f 

C 


(9) 


The radial eigenvalues ctj must "be determined from equation (8) and 
are a function of £/t) only. The methods used to determine the aj 
are presented in the appendix. 

Since the differential equation in pressure (eq. (3)) is linear, 
the pressure can he expressed as the summation of individual solutions. 
The pressure in the lined duct can then he expressed in a Fourier- 
Bessel series of solutions of the form of equation (5): 



= exp(iot) 



( 10 ) 


For an initial plane traveling pressure wave the pressure must he: 

P = exp(icub), (X = 0) (Oil) 

Combining equations (10) and (ll), the following equation for the 
constants Bj results : 

E b j j ° (“j t ) = 1 <12) 

The radial qharacteristic functions form an orthogonal set 
between r = 0 , and r = r w (ref. 5). The B^ can thus he deter- 
mined in equation (12), and the result is: 


B, = (13) 

■'"(M • 1 ] j o ( “j ) 

With equation (13) substituted in equation (10), the pressure any- 
where in the lined duct is determined. 

The complete solutions for the pressure and radial velocity also 
satisfy the radial impedance boundary conditions since: 




( 14 ) 


Axial Acoustic Energy. The axial intensity can he expressed 
as (ref. 6)1 
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(15) 


(16a) 


(16b) 

The subscripts jk signify the interaction of the j pressure 
mode with the k axial velocity mode. Equation (16a) must be used 
since: 

E jk / 0 (J / k) (17) 

In the evaluation of the integral of equation (16b) 

rJ o( a J t) J o( % t) dr 


+ C*)J 0 (af)J 0 (a k ) 

= (18) 

C£*(a| " a f } 

was used. 

Equation (18) was derived from the equation for the "Ortho- 
gonality of Characteristics Functions" given in reference 5. Equa- 
tion (8) and 

0o<O)0* = J 0 (ap (19) 

were also used in the derivation of equation (18). 

Equations (4), (5), (8), (13), (16b), and (18) yield: 



1 = 1/2 (Re (P#v x ) 

Axial acoustic energy can thusrbe determined from: 



4r w 

(0k+iT k )(C+^) e xp{-^[(a j +a k ) +i(T k - Tj)]x} 


^(“k - 1 J 


( 20 ) 


where 
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and Q* is determined from equation (21) by replacing k by j and 
taking 0 the complex conjugate. 

The damping and transmission parameters can be determined by sub- 
stituting equation (5) into equation (3). The real and imaginary 
parts of the resulting equation can be solved for cr. and r.. The 
result is : " 



( 22 ) 

where (+) is used for cr.* and (-) for t.. From equation (6) it is 
seen that Rj and 0j are the modulus and argumentof the complex 
eigenvalue aj. 

The sound power attenuation can be calculated from 


dB = 10 log]_Q 


^exit 

i 1 ' 

"entrance 


(23) 


where E eX i-fc a nd ^entrance are ca l cu l a ted from equation (16a) using 
equations (20), (21), and (22). Before equations (20), (21), and (22) 
can be evaluated, the complex eigenvalues aj must be calculated 
from the multiple solutions of equation (8). The method used to eval- 
uate the is outlined in the appendix. 

From equation (8) the radial eigenvalues aj are seen to be a 
function only of £/q. Equation (13) shows the J B., to be a function 
of £/q and a* or thus only of £/q. The modal 0 content of the 
plane pressure wave at the lined duct entrance (X = 0) is thus a 
function only of £/q according to equation (5). The damping and 
transmission parameters ctj and Tj are seen to be a function of 
£/q through Rj and 6y But, in addition, Oj and Tj contain 
q without ( explicitly. The exponential of equation X 20 ) contains 
the term expressible as: 



where L is the length of lined duct and D is the diameter. 

The variables on which the sound attenuation depends are thus 
£/q, t|, and L/D or alternately £, q, and L/D. 

Several pressure and velocity modes must be considered for soft 
walled ducts. The attenuation equation (eq. (23)) thus contains the 
Logarithm of the ratio of two complicated summations. Due to this 
complication no simplifying truncation of the theory yielding a simple 


8 


equation for attenuation has teen obtained. 
3. RESULTS AND DISCUSSION 


The numerical results which follow were calculated on an IBM 7094 
computer. The first ten radial modes were used. 

Attenuation Contours and Maximum Attenuation. The sound power 
attenuation depends upon four variables: the wall resistance and 

reactance, the duct diameter to wave length ratio (q), and the duct 
length to diameter ratio (L/D). Figure 1 presents constant attenu- 
ation contours as a function of wall resistance and reactance. The 
other variables are fixed at q = 1 and L/D = 3. Also plotted in 
figure 1 are the results of reference 1 which can be expressed as : 


dB = -17.36(L/D) ■ ' a * - (25) 

($ 2 + X 2 ) 

The following differences between the approximate theory of 
reference 1 and the more complete theory of this paper can be noted. 

For a given wall resistance maximum attenuation occurs at nega- 
tive wall reactance rather than at the wall tune point (X = 0). 

The present theory predicts an absolute maximum sound attenu- 
ation at fixed values of q and L/D. Decreasing resistance ulti- 
mately reduces attenuation for any value of reactance. The approxi- 
mate theory predicts local maxima, but for X = 0 the attenuation 
can.be increased indefinitely by reducing resistance. The present 
theory agrees with the approximate theory for large wall impedance 

(5/i). 

It should be noted that figure 1 is not a general attenuation 
plot which can be used for liner design. Figure 1 is valid only for 
a duct diameter to wavelength ratio ( q) of one and L/D = 3. For 
other q and L/D the impedance for maximum attenuation will move 
in the impedance plane, and the peak attenuation will change. 

Figure 2 presents the value of the absolute maximum sound atten- 
uation as a function of duct diameter in wavelengths ( q) . Several 
curves of constant L/D are shown. The maximum attenuation 
(-dB max ) has been normalized by dividing it by the value of L/D for 
each curve. The following should be noted from figure 2. For very 
small q (q «0.1) the maximum attenuation is weakly dependent upon 
q. Near ti = 1 the maximum attenuation falls rapidly with increas- 
ing q(a l/q 2 ). For large q maximum attenuation is approximately 
inversely proportional to q. 

For q > 1 the maximum attenuation is approximately proportional 
to L/D. An extrapolation of the larger L/D curves to q < 1 indi- 
cates that for small q the attenuation per unit L/D decreases 
with L/D. 

The locus of the maximum sound attenuation in the wall impedance 
plane is shown in figure 3. Several lines of constant q and L/D 
are presented. In general the locus for constant L/D shows that 
for increasing q (or f) the impedance for maximum attenuation moves 
to larger wall resistance and more negative reactance. The exception 
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is the inflection found in each constant L/D curve. The advent of 
this inflection occurs at roughly tj = L/D. 

The resistance of typical lining materials increases with fre- 
quency (see fig. 7). Unfortunately the reactance also increases with 
frequency. This trend of reactance with frequency is inverse to that 
required for maximum attenuation. Large sound attenuation can thus he 
obtained only over a limited frequency range for a given wall treat- 
ment. 


Sound Power Attenuation Contour Breadth in Wall Impedance Plane, 
Structural considerations may preclude the use of a liner material with 
optimum wall impedance. The question then arises. How far removed 
from the optimum wall impedance can the liner impedance be and still 
result in acceptable sound attenuation? Figures 4 and 5 are presented 
to partially answer this question. 

Figure 4(a) represents vertical cuts through the maximum attenu- 
ation points of several contour plots like figure 1. Figure 4(a) , for 
an L/D = 1, was constructed in the following way. For a fixed diam- 
eter to wavelength ratio (tj) and for the wall reactance of maximum 
attenuation, sound attenuations for several resistances on both sides 
of the maximum attenuation were calculated. Plots of attenuation 
versus resistance were made. These plots provided the resistances for 
say 5, 10, 20, etc. dB sound attenuations. These resistances were 
then plotted on figure 4(a) for the value of tj used. A range of tj 
was covered with the unique value of reactance for each rj which pro- 
vides maximum attenuation. Figure 4(b) was generated like figure 4(a) 
but with the roles of resistance and reactance interchanged. Fig- 
ure 4(b) thus represents horizontal cuts through the maximum attenu- 
ation points of several contour plots like figure 1. Figures 5(a) 
and (b) are analogous to 4(a) and (b) but with L/D = 5. 

The effect of L/D on the breadth of the sound attenuation con- 
tours is clear from figures 4 and 5, but numerical examples will 
further illustrate the effect. At the largest value of tj for which 
20 dB reduction can be obtained (tj = 1.13 for L/D = 1 and 
T) = 3.85 for L/D = 5) the breadth of the 10 dB contours were ob- 
served. For L/D = 1 the intervals are 0.2 < 0 < 1.6 and 

-1.41 < X < -0.15. For L/D = 5 they are 1.28 < $“< 8.9 and 
-6.5 < X < +0.5. The breadth increases (due to an Tncrease in L/D 
by a factor of 5) are thus 5.4 for resistance and 5.6 for reactance. 

The contour breadth is thus seen to increase approximately as L/D. 

The benefits possible from using lined partitions in the duct 
are now obvious. Increasing L/D by a factor of 5 reduces T] by 
a factor of 5 for the same L. Reducing t) increases the maximum 
possible sound attenuation (fig. 2). An increase in L/D broadens 
the impedance range over which a given sound attenuation can be 
obtained. This is caused by the increased liner area as well as the 
reduction in tj. 

Beaming of Pressure Wave at High Frequency . Radial sound pres- 
sure profiles for tj = 0.1, 1, and 10 for several values of L/D are 
shown in figure 6. The wall resistance and reactance were 1. 6 and 
zero for all curves of figure 6. 

At the lined duct entrance (L/D = 0) the pressure wave is a 
plane traveling wave of unit amplitude. This implies a uniform 
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amplitude and phase at any radius. At L/D = 3 the phase change 
from r = 0 to 0.9 r w is -4.6, -0.8, and +30.5 degrees for t) = 0.1, 

1, and 10 respectively. For large T] the pressure near the wall is 
seen to lead that at the duct axis. From these phase relationships 
and the radial pressure profiles of figure 6 the following is apparent. 
For t) = 0.1 and 1 the wave is attenuated and progresses down the duct 
as nearly a plane wave. For tj = 10 the pressure wave is beamed to- 
ward the duct axis. For L/D = 3 and 5 the pressure oscillations at 
the duct axis are more than twice that of the original plane wave 
(at L/D = 0). At the duct wall the pressure oscillations have "been 
reduced by a factor of about 10 (or 20 dB). The acoustic power has 
been redistributed toward the duct axis. At the wall where it must 
be absorbed, the sound intensity has been reduced. Thus, sound power 
attenuation is small as can be verified from the table of figure 6. 

Comparison of Theory and Experimental Results . In order to cal- 
culate sound attenuation in a duct, the resistance and reactance of 
the lining material must be known. The resistance and reactance fre- 
quency dependence (calculated by method of ref. 7) used in the follow- 
ing calculations is presented in figure 7. Several resistance curves 
are shown since liner resistance varies with sound pressure level in 
the duct. The following data were used in the impedance calculations: 

1. Steady flow velocity in duct, 350 fps. 

2. Liner open area ratio, 0.08 

3. Liner face plate thickness, 0.020 inches 

4. Liner hole diameter, 0.050 inches 

5. Liner back cavity depth, 1 inch 

Also included in figure 7 are two extra scales (r^ and 1^5) relating 
q(D/A) to frequency for L/D = 1 and 5 where L = 30 inches. 

Experimental sound power attenuation data was obtained with a 
lined inlet on a J-65 turbojet engine. The lined duct inside diam- 
eter (D) was 30 inches, and the lined length (L) was 30 inches (thus 
L/D = l) . The evaluation method of sound power levels is found in 
reference 9. The experimental sound power attenuation obtained with 
the lined inlet (properties as in figure 1 ) are shown by the data 
points in figure 8(b). 

Sound power attenuations (using the present theory) were calcu- 
lated for several sound pressure levels and a range of frequencies. 

The experimental apparatus dimensions were used. The results are 
shown in figure 8(a). The theoretical sound power attenuation is 
shown by the solid curve in figure 8(b). This curve was obtained by 
using the experimental sound pressure level at several frequencies 
and interpolating in figure 8(a). The agreement with experiment is 
seen to be excellent and within the experimental error (2dB). 

The dashed line of figure 8(b) shows the low frequency end of the 
attenuation curve predicted by the approximate theory of reference 1 
(see eq. (25) this paper). This curve continues rising with frequency 
and would predict about 18 dB sound power reduction at the linear tune 
point (4175 Hz). No reduction at this frequency was observed experi- 
mentally. 

The reason for the failure of this lined duct design to provide 
high, sound attenuation is seen in figures 4(a) and (b). The 130 dB 
level resistance and the reactance are plotted as dashed curves with 
frequencies noted on the curves. The liner resistance and reactance 
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diverge further from the region of high damping with increasing fre- 
quency. Above 1,000 hertz both resistance and reactance are in the 
low attenuation region. 

An example of how a geometry change can effect sound attenuation 
is shown in figure 9. The same liner material as in figure 7 was 
used. The lined duct length remains at 30 inches, but L/D was 
changed to 5 (say by lined partitions). The theoretical sound atten- 
uation was calculated for three sound pressure levels. High sound 
attenuation is seen to be possible over a wide frequency range. For 
comparison the 150 dB level sound attenuation curve of figure 8(a) has 
been replotted in figure 9. A change from L/D =1 to 5 is seen 
to increase the frequency for maximum attenuation from 1000 to 
2300 hertz and to increase the peak from 5. 3 to 51 dB. 

The 130 dB level resistance and the reactance from figure 7 are 
plotted on figure 5(a) and (t). A comparison with figure 4(a) and (b) 
shows how a change only in L/D from 1 to 5 has moved the high atten- 
uation impedance region. 

4. CONCLUDING REMARKS 


Comments on Duct Partitioning. Circular ducts have been par- 
tioned by the addition of radial lined sections (struts) and circum- 
ferential lined sections (rings) to increase the liner area to duct 
area ratio. With the proper selection of struts and rings the 
smaller resulting tubes can be made nearly square. Morse* s approxi- 
mate solution (ref. 1) for rectangular cross sections shows that 
equation (25) is valid for a square duct with very hard walls if D 
is replaced by a side of the square. Thus, in the limit of very hard 
walls the sound attenuation in circular and. square ducts is the same. 
For soft walls the analysis of this paper should give at least qual- 
itatively correct results for square ducts. 

Comments on Liner Resistance and Sound Power Level. As the in- 
itial plane pressure wave progresses down the duct it will be atten- 
uated and possibly beamed toward the duct axis (for D/A > l). The 
sound pressure level at the liner and thus the liner resistance 
(fig. 7) decreases down the duct. This effect could be handled by a 
new determination for the Bj of equation (10) at several locations 
along the duct. After the sound pressure level at the liner has drop- 
ped say 2 dB the characteristic solutions could be refitted to the re- 
maining pressure wave (not necessarily a plane wave). This new solu- 
tion could be used for a distance and then the refitting be repeated 
until the end of the duct is reached. 

Approximate solutions could also be divised. The extremes of the 
sound power level at the liner, for instance, could be used in an 
attempt to bracket the exact solution. 

The initial (duct entrance) sound pressure level was used in ob- 
taining the theoretical curve (solid curve) of figure 8(b). 

5. CONCLUSIONS 


A theory that describes the propagation and attenuation of an 
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initially plane traveling pressure Wave in a circular duct with an 
acoustically soft wall has been presented. The first ten character- 
istic function solutions were used in a Fourier-Bessel series to fit 
the plane traveling wave at the lined duct entrance. The results of 
the analysis may be summarized as follows. 

1. For a duct diameter to wavelength ratio greater than one the 
initial plane wave was found to beam toward the duct axis. This re- 
sults in a redistribution of the sound power within the duct. The 
sound pressure level at the wall liner, where acoustic energy must be 
absorbed, can be greatly reduced, while the sound power level is af- 
fected only slightly. 

2. The theoretical sound power attenuation was found to agree 
well with experimental data. The data was taken with a cylindrical 
duct, lined with a perforated plate material, mounted on the inlet of 
a J-65 engine. The experimental suppressor was designed on the basis 
of an approximate theory valid for hard walls. The results of the 
present analysis show the reason that the suppressor design failed 

to provide high noise reduction. 

3. For a given duct length to diameter ratio and duct diameter 
to wave length ratio an absolute maximum sound power attenuation 
was found in the liner impedance plane. This maximum was always 
found to be at negative wall reactance. With increasing duct diam- 
eter to wavelength ratio the maximum attenuation point was found to 
move to larger liner resistance and more negative liner reactance. 

The maximum sound power attenuation was found to be nearly propor- 
tional to L/D and to have a strong frequency dependence through 
the ratio of duct diameter to wavelength. 

4. The region in the wall impedance plane over which a given 
sound power attenuation can be obtained increases approximately in 
proportion to L/D. 

6. APPENDIX 


Calculation of the Bessel Function with Complex Argument . The 
calculation of the Bessel function with complex argument was done in 
a computer program subroutine. The following calculation methods 
were used. 

If R < 10 or if R < (1 + (20/jt))n the infinite series defin- 
ition (ref. 5) of the Bessel function was used as: 


or 
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J n (Re 10 ) 


(-D 




k/Re 


i0\ 


2k+n 


2 ) 

kJ (k + n).' 


k=0 
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} (26) 
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2k+n 


J n (R ei9 ) 


T^o 


J 


(-D k (-\ 

;(k + n j ?" E os ( 2k + n ) 0 + 1 sin(2k + n)e) 
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If R > 10 and R > (1 + (20/jr))n, the following asymptotic 
expansion (ref. 5) was used: 


J n (Re ie ) ~ ^-^16 {u n (Re i0 )c'os[G n (Re i0 3 - V n (Re i0 )sin[G n (Re i0 )]J 


where 


G n (Re 10 ) = Re i0 ~ j ~ ~ 


(27) 

(28) 


U n (Re 10 ) = 1 + 


00 

E 

k=l 




(- 1 ) 

(2k)»(8Re i0 ) :iK 


and 


V n (Re i0 ) = 


S. = (4n^ 
2 


00 

£ 

k=l 


/ x k-l 
(-1) S 


4k-3 


> 


(2k - l) , (8Re ie ) 


2FT 


l 2 )(4n 2 - 3 2 ) ( 4n 2 - 5 2 ) . . . (4n 2 - j 2 ) 


(29) 


The series of equation (29) were terminated in two ways: when 

the modulus of the k th term of U n (Re ie ) exceeded the modulus of 
the k^ h term of V n (Re i0 ), or when the modulus of the k^* 1 term 
of either was less than lO - ^. 

Rapid convergence occurs in equation (26) for R « n and in 
equation (27) when R » n. Both equations (26) and (27) were tested 
over wide ranges of R, 9, and n. For large n the dividing line 
between the two methods is R = (1 + (29/rt))n. If R < 10, equa- 
tion (26) converges more rapidly for any n. 

Double precision arithmetic was used in all cases. 

Calculation of Complex Eigenvalues. The complex eigenvalues 
must be determined from the solution of equation (8) (derived in 
the text and repeated here). 


Z _ £ _ iitJ 0 (cc,j) 
Po C t 1 n ” O'jJi(aj) 


where again: 


aj - R^ 0 3 


( 8 ) 


( 6 ) 


The method of reference 4 was used in the solution of equation (8) 
for large wall impedance. 
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Let 



Using equation (6), substituting equation (30) into equation (8), and 
then differentiating equation (8) with respect to e, there results: 

(F + € 2 ) || = F (31) 

Equation (3l) is a first order nonlinear differential equation which 
can he solved for F in an infinite series solution in e. 

The boundary conditions which must he applied are: 
at e = 0 

= 0 (all J) 

and Rj is the root of J-^(Rj) = 0 

or Rj = 0, 3. 8317, 7.0156, etc. for j = 1, 2, 3, etc. 
d For small wall impedance let 

! = | (33) 



Substitution of equation (33) into equation (31) yields: 

(1 + | 2 F) M + F = 0 (34) 

d| 


The boundary conditions are: 


at |=0 0j = 0 
and Ri is the 

O 


(all j) 


root of J q (R.) = 0 


(35) 


or R^ = 2.4048, 5.5201, 8.6437, etc. for j = 1, 2, 3, etc. 

For every mode there is a region in the impedance plane where 
the series solutions of both equation (31) and equation (34) diverge 
(at least on the computer using double precision arithmetic). This 
region is around the branch line. The branch line is defined here 
as the curve in the liner impedance plane ($/t], X/t]) on which two 
solutions to equation (8) have the same R but different 0(Rj = % 

but j 9 k ). The critical values of R. on the branch line are 

3.1962, 6.3064, ... or: d 


FjjCrit *** J 11 (36) 

In this critical impedance region the method of reference 8, for 
solving nonlinear equations, was applied directly to equation (8). 
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Figure 4(a). - Sound power attenuation contours at reactance of maximum 
attenuation (fig. 3) and L/D = 1. 



Figure 4(b). - Sound power attenuation contours at resistance ot maximum attenuation (fig. 3) and L/D * 1. 
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Figure 5(a). - Sound power attenuation contours at reac- 
tance of maximum attenuation (fig. 3) and L/D = 5. 



Figure 5(b). - Sound power attenuation contours at resistance of maximum attenuation (fig. 3) and L/D - 5. 







